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Abstract 

We consider the Darboux transformation of the Green functions of the 
regular boundary problem of the one-dimensional stationary Dirac equation. 
We obtained the Green functions of the transformed Dirac equation with the 
initial regular boundary conditions. We also construct the formula for the 
unabridged trace of the difference of the transformed and the initial Green 
functions of the regular boundary problem of the one-dimensional stationary 
Dirac equation. We illustrate our findings by the consideration of the Dar- 
boux transformation for the Green function of the free particle Dirac equation 
on an interval. 
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1 Introduction 

There has been great interest in using there Darboux transformation [H [2] for the 
analysis of physical systems [3], HI |5] and for finding new solvable systems [6], [71 
[H]. It has been shown that the transformation method is useful in finding soliton 
solutions of the integrable systems [HI [ID] and constructing supersymmetric quantum 



1 



mechanical systems [HI [T2] o It is well known that SUSY QM is basically equivalent 
to the Darboux transformation and the factorization properties of the Schrodinger 
equation p ^lT^flS] . The Darboux transformation of the one-dimensional stationary 
Dirac equation is equivalent to the underlying quadratic supersymmetry and the 
factorization properties [T^ \T7\ flSj. 

Despite of the growing number of papers in this field many questions still remain 
open and require further study. In particular, the author are aware of only several 
papers [IHl EHl El] devoted to the SUSY transformations at the level of the Green 
functions. In [T9] was obtained the integral relation between the Green functions 
for two SUSY partner Hamiltonians of the one-dimensional Schrodinger equation 
with discrete spectra. In [201 EI] the integral relation between the Green functions 
corresponding two SUSY partner Hamiltonians of the one-dimensional Schrodinger 
equation is generated to the case of continuous spectrum. The exact Green function 
of the time-dependent Schrodinger equation was obtained in [22l [23] . 

An interesting open problem is to find the analogous results for the Dirac equa- 
tion. In this paper, we construct the Darboux transformation of the Green function 
for the regular boundary problem of the one- dimensional stationary Dirac equation 
with a generalized form of the potential and obtain formulas for the unabridged 
trace of the difference of the modified and initial Green functions. The rest of the 
paper is organized as follows. In Section 2 we construct a Green function for the 
initial regular boundary problem of the one-dimensional Dirac equation with the 
generalized form of the potential. In Section 3 we consider the Green function for 
the transformed regular boundary problem of the one-dimensional Dirac equation. 
We construct the Darboux transformation of the Green function. In Section 4 we 
obtain formulas for the unabridged trace of the difference of the modified and the 
initial Green functions and consider the spectral representation of the correspond- 
ing unabridged trace. In Conclusion the summarize our results and speculate about 
some perspective. 

2 Green function of the one-dimensional Dirac 
equation 

The Green function of the one-dimensional Dirac equation 

(hoix) - E)i;{x) = (1) 

with the Dirac Hamiltonian of the form ho = ia2dx+V{x) ^24] is needed for obtaining 
the solution ^ 

<!>{x,E)= [ Go{x,y,E)F{y)dy (2) 

J a 

of the following inhomogeneous equation 

{ho{x)-E)<l>{x,E) = F{x), (3) 

"'^ Recently a special issue of Journal of Physics, A 34 was devoted to research work in super- 
symmetric quantum mechanics (SUSY QM). 



2 



where a and b are the endpoints of an closeg interval. 

Earher this method allowed the authors of the paper [25] to obtain the Green 
function of the one- dimensional Dirac equation with the potential 



Un 



p{x) 
q{x) 



In the present paper, we will obtain the Green function of the Dirac equation 
with the generalized form [2S] of the potential 

Vo{x) = u{x)I + {m + S{x))cx3 + q{x)cxi, (4) 

where uj{x), S{x) and q{x) are real functions of x, m is the mass of a particle, Ui, 
(Ts are usual Pauli matrices. The generalized potential is the self-adjoint potential. 

We consider the following boundary conditions for the components $i (x, E) , 
$2(2;, of the solution (ED to ([3]): 

$i(a, E) sin(a) + $2(0, E) cos{a) = 0, (5) 



$1(6, E) sm{(3) + $2(6, E) cos(/3) = 0. (6) 

We suppose that the solutions iplx), ip{x) of the Dirac equation with the gener- 
alized form of the potential obey the following conditions: 

/ f cos a \ , /, f COS 3 \ 

^^«'^)=(-sina J' ^^^'^)=(-sin/3 J" 

Let us construct the matrix 

r (r V F\- ^ \ ^lj{x,E)ip^{y,E), y<x, . . 

°^ I x<y, 

where W{E) = W{ip{x, E),iIj{x, E)} = const is the Wronskian of the two functions 
ipi^x) and (p{x). 

One can reality examine that the vector-function ^{x,E) is the solution to the 
inhomogeneous equation and therefore the matrix ([H]) is the Green function of 
the regular boundary problem of the one- dimensional stationary Dirac equation with 
the generalized form of the potential (jll). 

T 1 • / r-i\ Vli^jE) 

Let us construct the spmor yix, E) = \ ; „; 

y{x,E)= [ Go{x,y,E)F{y)dy. (9) 
Jb 

In detail, Eq. ([H]) becomes 



Lpi{x,E)ipi{y,E) Lpi{x,E)ip2iy,E) 

r „ v\- " ) \ Mx,E)i;i{y,E) ^p2{x, E)ip2{y, E) J ' , . 

^o[x,y,i^) ^^^^<| ^ ^^(^^^)^^(^^^) cp2iy,E)Mx,E)'^ ^'^^ 

Vi{y,E)ilj2{x,E) ip2{y,E)ilj2{x,E) 



y < X. 
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We write down a detailed form of the components of the spinor 

r2 

1 ( Mx,E)My,E)F,{y)+ipiix,E)My,E)F2{y), x<y 



Go{x,y,E)F{y)=[ ), (11) 



W{E) \ Vi{y, E)iIj^{x, E)Fi{y) + ^^{y, E)F2{y), y < x, 

P 1 / M^,E)My,E)F,{y)+ip2ix,E)My,E)F2{y), x<y 

' W{E) 1 ipiiy, E)iIj2{x, E)F,{y) + ^^{y, ^)^2(a;, E)F2{y), y<x. 

From these equations, taking into consideration (jlip and 

E)FM + My. E)F2iy) = E)Fiy), 

E)Fr{y) + ^y, E)F2{y) = '/{y, E)F{y), 

we find the components of the spinor ^ yi{x,E), y2{x,E). 

We consider the components of the spinor ([9]) in the case y <x, then y G [a,x], 
after that in the case y > x, when y G (x, b], and as a result, we find the components 
of the spinor ([9]) yi{x,E), y2{x,E) in the general case: 

yi{x,E) = ,J^^{Mx,E)J^ ip^{y,E)F{y)dy+ 

+ ipi{x,E) [ ^'^iy,E)F{y)dy}, (12) 



Z/2(x, E) = :^^{^2(x, i^;) ^ </^^(i/, i?)F(|/)f/|/+ 

+ yP2(x,i?) /" ^''{y,E)F{y)dy}. (13) 

Now we check that the vector-function y{x, E) is the solution of the inhomoge- 
neous equation ([3]) 

dM^, E)-{E- S{x) - m - uj{x))yi{x, E) + q{x)y2{x, E) = F,{x), (14) 

d,yi{x, E) + {E + S{x) +m- u;{x))y2{x, E) - q{x)yi{x, E) = -F2{x), (15) 

where Fi{y), F2{y) are the components of the spinor F{y). 

Let us show that equality ( |T5l) is true. For that we differentiate yi{x,E): 

d,yi{x,E) = ^J^{d,Mx,E) ip^{y,E)F{y)dy 
+ d,Mx,E) / ij''{y,E)F{y)dy} 

J X 

- ifiix, E){Mx, E)Fi{x) + V'2(x, E)F2{x)} 

+ V'i(x, E){^^{x, E)F,{x) + M^, E)F2ix)}. (16) 
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Since the spinors iplx, E) and (p{x, E) are the solutions to the Dirac system of 
equations, the following equalities are correct: 

d^%l)i{x,E) = q{x)ilJi{x,E) - {S{x) +m-uj{x) + E)'^2{x,E), (17) 



d^ipi{x,E) = q{x)(pi{x,E) - {S{x) +m - uj{x) + E)ip2{x,E). 
Substituting equahties (fTTl) . (|T8l) into (fT6l) . we obtain 

d,yi{x,E) = ,^^[q{x){Mx,E) ip^{y, E)F{y)dy 
+ Mx,E) [ tP^iy,E)F{y)dy} 



- (Six) +m- uix) + E){^/J2ix, E) / (y, E)F{y)dy 



a 



(19) 



+ ip2{x,E) [ tP^{y,E)F{y)dy}]-F2{x). 

From ([12]), dni) for yiix,E), y2{x,E) we get 

9,1/1 (x, E) = q{x)yi{x, E) - {S{x) + m - u{x) + E}y2{x) - F2{x). (20) 

Therefore, ( JTSl) is correct. Similarly, equality ( fMl) can be proven. The spinor- 
function $(x, i?) = y{x,E) is the solution to ([3]) with the generalized form of the 
potential. 

Now we demonstrate realization of the regular boundary conditions for 
spinor-f unction $(x, E) = y{x, E) 

yi{a, E) sin(Q;) + ^2(0, E) cos(a) = 0, (21) 

yi(6, E) sin(/5) + y2{h, E) cos(/5) = 0, (22) 

where yi{x,E), y2{x,E) are the components of the spinor y{x,E). We compute 
(fT2l). fffsll when x = a, x = b: 

yi{a,E) = ^^^{i;^{a,E) ^'^ {y, E)F{y)dy 

+ ^i{a,E) I 'iP^{y,E)F{y)dy}, (23) 



y2{a,E) = ^^^{,l,2{a,E) I y,^{y, E)F{y)dy 



a 



+ ^2{a,E) j\^{y,E)F{y)dy], (24) 

J a 
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yi{b,E) = ,^^{Mb,E) I ^^{y,E)F{y)dy 



+ ^,{b,E) f\^{y,E)F{y)dy}, 
Jb 



(25) 



y2{h,E) = ^^^{i;,{h,E) j\^{y,E)F{y)dy 



b 

T 



+ Mb,E) / riy,E)F{y)dy}. (26) 
Jb 

Since f{y)dy = and f{y)dy = 0, where f{y) is the arbitrary function of y, 
the equahties ([23]), (El, (fe]), (EHj) take the following forms: 

yiia, E) = E) ^ ^^(y, E)F{y)dy}, (27) 

|/2(a,i?) = ,^^{^,(a,E) iy^{y,E)F{y)dy}, (28) 



^) = ]^^iMb, E) I v,^(y, E)F{y)dy}, (29) 



l/2(&, i^^) = ^^^{Ub, E) ^ y.^(y, i5;)F(y)rfi/}. (30) 

Now, taking into account (I2ID> (EHD, (EHD, dSDl), we check the conditions ^ 
yi{a, E) sin(a;) + 2/2(0., E) cos(a) 



1 



W{E) 



b 



{<^i(a,E)sin(a) + (^2(a,^)cos(a)} / (y, E)F{y)dy, (31) 



yi{b,E) sin(/3)+ 2/2(6,^) cos(/3) 
1 



{^1(6, E) sin («) + V^2(&, E) cos (a)} / (/.^(y, E)F{y)dy. (32) 



Due to the boundary conditions ([7]) for the functions 'ip{x, E), ip{x, E) the relations 
are (|2T1) . ( !22l) are valid. Hence, the matrix Go{x,y, E) (El), is the Green function 
for the regular boundary problem of the one-dimensional stationary Dirac equation 
with the generalized form of the potential. 

The matrix ([H]) can be written in the following form: 

Go{x, y, E) = (^(x)<^^(y)e(x -y) + ^{xW {y)Q{y - x)) / {W {^{x),^{x)}), (33) 
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where <d{x — y), Q{y — x) are the Heaviside step functions. Thus, we have con- 
structed the Green function of the initial problem. 



Example 1 Let us calculate the Green function of the initial problem for the Dirac 
([T]) free particle equation with the = ia2dx + Vq{x), Vq = ma^. 

We first chose the two linearly independent solutions ip, ip for the free particle 
case: 

, _ / cos (A;x) \ _ / cosA;(x — 1) \ , 

^ ~ \ksm{kx)/{E + m) i' ~ \ ksiiik{x - 1)/{E + m) J' 



k'^ = E^-m^, W{(fi,ip} = -k sin {k)/{E + m). (35) 

Next we apply the boundary conditions to functions ip, y9 for the case a = 0, 
6=1: 

^ \x=i = ± J ^ , \x=o= ± J ^ . (36) 

From ( l34l) it follows that k = im, n & Z . 

Also, we find the final expression for the Green function of the initial regular 
boundary problem ([5]), with a = P = ivrn: 

r I p^ A^Q{x-y) + B^Q{y-x) 



/ COS {kx) cos {ky — k) k cos {kx) sin {ky — k)/{E + m) 

k sin {kx) cos {ky — k)/{E + m) k"^ sin {kx) sin {ky — k)/{E + rnf' 

= {AY{x ^ y). 

3 Darboux transformed Green function 

Before constructing Darboux transformed Green function, let us shortly consider 
the Darboux transformation method for the Dirac equation. 

The so-called Darbuox transformation operator L and the potential Vi of the 
transformed Dirac Hamiltonian hi have respectively the following form: 

L = d^-UxU'\ (38) 
Vi = Vo + [ta2,u,u~']. (39) 

The so-called transformation function u is 2 x 2 matrix consists from two solutions 
u = {ui,U2) of the initial Dirac equation with E = Ai,A2, where Ai, A2 are the 
neighboring energy levels. 
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If the function if) is solution of initial equation and ip 7^ mi, M2, then the function 
if) = Lip is solution of the transformed Dirac equation 

{hi{x) - E)tP{x) = 0, hi=ia2 + Vi. (40) 

The Green function of the equation (HOl) with the regular boundary conditions 
for the components $i(a;,-E), $2(2^, -E) 

l>i(a, E) sin(a) + $2(0, E) cos(a) = 0, (41) 



$1(6, E) sin(/3) + ^2{h, E) cos(/3) = (42) 
of the spinor ^{x,E) = j^Gi{x,y,E)F{y)dy one can represent in the form: 
Gi(a;, E) = {^{x)(fF{y)Q{x - y) + 0{xW {y)Q{y - x)) / {W {^{x),^{x)}). (43) 

(i) If (5 = a, /5 = /3, we have the Green function of the Dirac equation with the 
transformed potential and the initial regular boundary conditions. 

(ii) If a 7^ a, /5 7^ /5, we have the Green function of the Dirac equation with the 
transformed potential and the modified regular boundary conditions. 

Further we consider only the case when a = a, (3 = (3. 

Now let us consider the construction of the Darboux transformed Green function 
on the basis of the concrete example. 

Example 2 Let us first chose the transformation functions as follows: 

_ / 1 cos {kix) \ 

" ~ VO A;isin(A;ix)/(A + m) J' ^ > 

Next, obtain from the free particle case the transformed potential 

Vi = — AfTs + fci cot (A;ix)(Ji, (45) 

where 



Ai = m, A2 = A = ±Vm2 + 7r2. (46) 

It can be checked that only these two transformation matrices (differing from each 
other only by sign of A) provide the conservation of the initial regular boundary 
condition. 

Now the Darboux transformed functions tp, (p one can represent in the form: 

~ ^ {\-E)smkx \ 

E + m \ k cos kx — ki cot kix sin /cx y ' 

^ ( (A-E)sinfc(x-l) \ 

E + m \ k cos {kx — k) — ki cot (kix) sin {kx — k) J ' ^ ' 
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Figure 1: The function f{x) = ki cot (fcix), ki = n 



Finally, the Green functions of the Dirac equation with the transformed potential 
)]) may be given by the formula: 

. A'"9(x-^+i^W9(,_,) ^ ^^^^ 

W{if,i^} = {E - X){E ~m)W{ip,ilj}, (50) 



where 

- ( S 1 1 • - '■^'")^<^ " <^^' 

aS\^ = (A-E)2sin(A;x)sin(%- A;), 

^21'' = (A — £')(/(; cos (/cx) — /ci cot (/cix) sin (fca;)) sin (% — A;)), 

= (A — i?) sin (A;x)(A;cos (% — A;) — /ci cot (/ci?/) sin (A;?/ — /c)), 
^22^ — cos (A;x) — ki cot (A;ix) sin (kx)) 

X (A; cos (A;?/ — k) — ki cot (A;i?/) sin (A;?/ — A;)). 

Thus, we have constructed the Green function of the transformed Dirac equation 
with the potential and the initial regular boundary conditions. 



4 Unabridged trace for the difference of the 
modified and the initial Green functions 

In this section we calculate the unabridged trace of the difference of the transformed 
and the initial Green functions. 
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Firstly, let us consider the unabridged trace of the transformed Green function 
tr Gi{x,y, E)\x=ydy. It can be represented in the form: 

tr Gi{x,y,E)\^=ydy = ° • (52) 

Ja Wi^^ij} 

It is known that the action of the transformation operator and the conjugate trans- 
formation operator on the spinors may be written in the following way 



= u-^(u-V), (53) 
dx 

L+xP = -(n+)-i-^(M+^). (54) 



Also, accounting these properties, we can rewrite flozl) as: 



tr J Gi{x,y,E)\^,=ydy = tr J —{u ^lj)(p u\x=ydy ^ Wiifip}. (55) 

Next let us consider the unabridged trace of the difference of the transformed 
and the initial Green functions. It is obvious [23] that the spectral representation 
for the initial and the transformed Green functions correspondingly look like as 

Go{x, y,E) = 2_^ ^ _^ — , (56) 

no 

Gi{x,y,E) = E -E ' ^ ^ 

m " 

where the functions 

= {{En, - Ai)(K, - X2))-'L^n,L, (58) 



form an ortho normal set. The completeness of (1581) will be shown in examples 3 and 
4. 

The construction 

tr [ {Giix,y,E)-Goix,y,E))\,=ydy (59) 

J a 

is called the unabridged trace of the difference of the transformed and the initial 
Green functions. The spectral representation of fl59|) is as follows: 

tr j iG,{x,y,E)-Goix,y,E))U=ydy = J2 \e ~^ E - 

ni no 

where En^ are discrete eigenvalues of Hq and hi respectively. 
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If the set (1581) is complete then the formula (!59l) one can represent in the form: 

tr ^ y, E) - Go(a;, E))l=ydy = + (61) 

Now let us show an another derivation of the trace formula for the Green func- 
tions difference. Let us integrate (1331) by parts and apply the trace property 

trAB = trBA (62) 

to obtain the relation 

tr \ G,{x,y,E)\,=ydy= ''l^Z- ■ (63) 



Due to ( |50l) and the factorization property from [16] we obtain: 

/ trGi{x,y,E)\,=ydy = JJ^ '° + / trG^{x,y,E%=ydy. 

Ja Wiipip} Ja 

Similarly, we would like to write 



(64) 



/ trGi{x,y,E)\^=ydy = J" + / trGo{x,y, E)\^=ydy. (65) 

Ja W{(ptp} Ja 

Finally, from flU^ . (1651) we find: 



tr [ {G^{x,y,E)-Go{x,y,E))U=ydy = }_ tr{ij^^) t (66) 

J a W{(p,1p} 

tr f iG^{x,y,E)-Goix,y,E))U=ydy = ] tr(^vp^) |^ (67) 

J a W{ip,1p} 

Here a, 6 are limits of the integration. As shown in Appendix, right parts of (1661) . 
( 1671) are equal. 

Example 3 In this example we calculate (l59ll with the help of ( l67j) . 

More precisely, take into account (!35ll . (HTll . (|50l) . we calculate tr(-?/'y9"^)/iy{(^, ■?/'} 

a 

where a = 0, 6 = 1: 
1 



tr{'ijj(p ) 



(A — ii^) sin(A;x) cos(A;x — k) + 



{E - Tn){E - X)sm{k) j \ j \ ^ E + m 



cos(fcx) sin(A;a; — k) — cot(/cix) sin(fca;) sin(A;a; — A;) 

E + m 







Using the L'Hospital's rule for evaluating of indeterminations, we obtain 



tr{%pip^ 
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E~m E~\ 
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Due to (1671) we get: 

tr [ {Gi{x,y,E)-Go{x,y,E))\,=ydy = -^ + -^, (69) 

Jo — Ai tj — \2 

where Ai = m, A2 = A. This result indicate that the spectrum of the transformed 
Dirac Hamiltonian differs from the spectrum of the initial Dirac operator. The 
spectral lines E = X,m disappear from the transformed spectrum and the set of 
functions </)„j is complete. 



Example 4 In this example we calculate (|59|) in a direct way. 
First, we consider 

^ / E mcos(2A;x — A;) , , 

trGoix, x,E) = -- cot A; , ^ , ^ 70 



Next, we integrate the expression ( 1701) 

/"I ^ 

t''" / Go(x, X, E)dx = —— cot k — 
Jo k k^ 

k = Trn, cot{k) - k~^ = i:'^^^2k/{E^ - El), 

tr [ Go{x,x,E)dx = -^n=i{E + E^)-'^ -J:^^^{E - En)~^ - {E -m)-\ 
Jo 

Similarly, we calculate 
»i 

v'oo / 771 771 \ — 1 I / 771 \\— 1 



(71) 



(72) 



tr / Giix,x,E) = -J:^=i{E + Er.)-'-J:Z,{E-E^)-' + {E-X)-' (73) 
Jo 

and obtain that 

f 11 

tr Gi{x,x,E)~Go{x,x,E)dx = - + - . (74) 

J E — A E — m 

Thus, we can conclude that the set of functions (pm is complete. 



5 Conclusion 

In this paper, we have studied the Darboux transformation of the Green functions 
of the regular boundary problem corresponding to the initial and the transformed 
potentials of the one-dimensional Dirac equation for the case of the Dirac Hamilto- 
nians with discrete spectrum. The main results of the paper are the construction of 
the Darboux transformed Green function with initial regular boundary conditions 
and the formulae for an unabridged trace ( l66i) . ( 1671) . For the checking of formulae 
(1661) and (1671) we consider the unabridged trace of difference between transformed 
and initial Green functions by usual possible and by formula (lUUj) . both results are 
equal. The all results of this paper are studied only for discrete spectrum and regu- 
lar boundary problem. An interesting question for the future is generation of these 
results to the case of continuous spectrum of the Dirac equation and the case of prob- 
lem on the real line and a half-line. We believe that these problem will investigate 
in a separate publication. 
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Appendix. Equality of traces 

In this Appendix, taking into account ( IMl) . ( l35i) . (H71) . (HHI) . we show the equivalence 
of expressions for right parts of (1661) . ([6] 



tr{i'(fy')t = {i'i^i + ^2^2)t=i^i{h)^i{h) 

+i)2{h)^2{h) - i)i{a)(pi{a) - i)2{a)(p2{a), (75) 

+^2{b)'f2{b) - i'iia)ipi{a) - 4j2{a)'f2{a), (76) 

where a = 0, 6 = 1. 

In the explicit form (1751) . ( 1761) look like as follows: 

tr{ipLp'^) Iq = {(A — E) cos{kx) sin(A;x — k) + {E — m) sin(A;x) cos(A;x — k) 
kki 

— — cot(A;ix) sin(A;x) sin(A;x — A;)} li, (77) 

E + m 

tri^xpLf'^) Iq = {(A — E) sin(A;x) cos(A;x — k) + (E — m) cos{kx) sin(A;a; — k) 
kki 

cot(A;ix) sin(A;x) sin(A;a; — k)} li . (78) 

E + m 

Since 

tr(^^^) 1^ -tr(^^^) \l = ((A - E) sin(-A;) + {E - m) sm{k)) \l= (79) 
we finally obtain 

trii:^^) \l= tr{^^^) \l . (80) 
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